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INTRODUCTION 


Common  testing  procedures  for  the  laboratory  measurement  of  material 
properties  for  use  in  ground  motion  calculations  have  generally  consisted 
of  standard  hydrostatic,  uniaxial -strain  and  triaxial  tests.  It  has  recently 
been  recognized  that  these  paths  are  not  necessarily  the  ones  that  are 
followed  in  actual  field  applications,  i.e.,  conventional  and  nuclear 
explosions  in  the  earth.  Since  difficulty  is  often  experienced  in  developing 
accurate  constitutive  models  that  are  valid  for  a  wide  range  of  loading 
conditions,  it  seems  important  to  follow,  as  closely  as  possible,  the  stress 
paths  (or  strain  paths)  that  are  experienced  by  material  elements  in  actual 
field  conditions.  Furthermore,  since  measurement  techniques  do  not  yet 
allow  the  field  determination  of  these  stress  paths  (or  strain  paths),  one 
must  rely  on  numerical  calculations  and  an  initial  best  estimation  of  the 
material  constitutive  properties.  In  this  report  we  present  the  results 
of  one-dimensional  numerical  finite-difference  calculations  for  cylindrical 
and  spherical  wave  propagation,  which  define  the  stress  and  strain  paths 
followed  by  material  elements  at  varying  distances  from  cylindrical  and 
spherical  explosive  sources  in  the  earth.  The  purpose  of  these  calculations 
is  to  define  laboratory  tests  best  suited  for  the  definition  of  materia'l 
constitutive  behavior  in  the  analysis  of  CIST  (Cylindrical  In  Situ  Tests) 
and  other  subsurface  explosive  events.  On  the  basis  of  these  calculational 
results,  static  laboratory  tests  are  conducted  which  represent  strain 
paths  experienced  by  material  elements  in  the  vicinity  of  cylindrical  and 
spherical  explosions  in  an  infinite  medium.  The  material  tested  in  the 
experimental  program  is  Ka.yenta  sandstone. 


STRESS  PATH  DETERMINATION  FROM  FINITE-DIFFERENCE  SOLUTIONS 

The  quantities  which  are  obtained  from  the  finite-difference  solution 
are  ck  and  as  functions  of  time  at  various  distances  from  the  explosive 
source.  For  purposes  of  definite  laboratory  tests,  it  is  useful  to  express 
the  output  of  these  calculations  in  terms  of  the  load  L  =  o  -  p  and  p 
in  the  triaxial  test  configuration.  Here  aa  is  the  axial  stress  and  pc  is  the 
confining  fluid  pressure.  It  is  also  more  convenient  to  deal  with  axial 
and  transverse  strain  components  (e_  and  e+)  in  the  triaxial  test  rather  than 

a  t 

e.j  defined  in  the  finite-difference  solution.  In  the  case  of  spherical 
flow,  one  would  simply  make  the  identification  that  L  =  aj  -  o3, 

Pc  =  o 3,  ca  =  ej,  and  =  e3.  For  cylindrical  flow  the  identification 
is  slightly  more  complicated. 

In  general,  let  us  assume  that  we  have  values  of  stress  and  strain 
invariants  defined  by 


t(t)  2  [ (a i"02 ) 2  +  (°2-a  3 ) 2  +  (oa-Oi)2]1^2//^ 

(1) 

p(t)  =  (aj  +  o2  +  o2)/3  , 

(2) 

ev(t)  -  cj  +  c2  +  £3 

(3) 

ed(t)  £  [(ei-c2)2  +  ( e2~ e 3 ) 2  +  (e3-Gi)2]1/,2/>^  > 

(4) 

as  functions  of  time  at  a  fixed  spatial  position  as  provided  by  the  finite- 
difference  calculation.  If  the  material  constitutive  behavior  involves 
only  first  and  second  invariants  of  the  stress  and  strain  tensors,  the 
quantities  defined  by  Eqs.  (1)  -  (4)  can  also  be  written  in  the  following 
terms  for  the  purpose  of  defining  laboratory  test  paths: 
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e  (t) 
V 


ed(t) 


(oa  -  Pc)/^3 


(°a  +  2pc)/3 


p  + 

a  ct 


(£a  ‘  et)/,/3 


and  hence  laboratory  stress  and  strain  paths  become  in  parametric  form 
(t  as  the  parameter): 


L  =  /3  x(t) 


Pc  »  P(t)  -  t( t)/i/3 


ea  =  S(t)/3  +  2ed(t)//3 


et  -  ev(t)/3  -  E(J(t)/^ 


Calculational  Results 

Stress  (and  strain)  paths  for  cylindrical  and  spherical  wave  propagation 
have  been  calculated  with  the  use  of  elastic-plastic  constitutive  descriptions 
presented  in  the  Appendix.  The  material  parameters  are  chosen  to  be 
representative  of  Mixed  Company  (Kayenta)  sandstone.  In  all  cases  a  radial 
stress  given  by 


0-  =  P«e 
r  o 


is  applied  at  the  interior  cavity  surface  of  radius  RQ  =  1  m.  The  peak 
radial  stress,  pQ,  is  taken  to  be  10  kbar  and  the  decay  constant,  1/a,  takes 
on  values  of  0.1  msec,  1.0  msec  and  10  msec.  All  results  are  presented  in 


terms  of  c  vs.  l.  (axial  strain  vs.  transverse  strain)  and  L/u  vs.  p  /K 

3  U  C 

(load/shear-modulus  vs.  confining-fluid-pressure/bulk-modulus) ,  i  .e. ,  the 
quantities  related  directly  to  static  triaxial  laboratory  tests. 

Figures  la,  lb  and  1c  show  stress  and  strain  paths  at  various  distances 
from  a  cylindrical  explosion.  At  the  radial  position  R  =  2Rq  the  stress 
path  intersects  the  failure  surface  during  loading  and  remains  in  contact 
during  unloading.  The  corresponding  strain  path  initially  approximates 
conditions  of  uniaxial  strain,  but  exhibits  considerable  transverse  strain 
during  the  latter  stages  of  deformation.  At  R  =  3Rq  it  can  be  seen  that 
the  strain  path  is  approximated  by  loading  in  uniaxial  strain  followed  by 
unloading  at  constant  axial  strain,  while  at  R  =  5Ro  the  axial  strain  is 
seen  to  decrease  during  unloading.  Of  course,  at  much  greater  distances 
from  the  explosive  source  plane-wave  conditions  are  achieved,  and  the  load- 
unload  path  remains  on  the  et  =  0  axis. 

Figures  2a  -  2e  show  similar  behavior  for  spherical  wave  propagation. 
Figures  1  and  2  give  an  indication  of  how  strain  and  stress  paths  depend 
on  distance  from  the  source.  Another  important  consideration  is  that  of 
pulse  shape  or  pulse  duration.  This  is  controlled  by  the  parameter  a  in 
Eq.  (13).  A  number  of  calculations  were  performed  for  cylindrical  geometry 
with  1/a  =  0.1  msec,  1.0  msec  and  10  msec.  The  peak  radial  stress  pQ  remains 
the  same  in  all  calculations  (pQ  =  10  kbar).  The  resulting  stress  and  strain 
paths  are  shown  in  Figure  3  at  radial  positions  1.5R  ,  2R  ,  3R  ,  4R  and 
5Rq.  One  sees  immediately  that  not  only  does  the  position  have  influence 
on  stress  and  strain  paths,  but  also  that  pulse  duration  has  a  significant 
effect.  It  will  therefore  be  important  to  represent,  as  accurately  as 
possible,  the  time  history  of  the  cavity  stress  due  to  the  explosive  source. 
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et  and  pc/K  (x  10  ) 

Figure  la.  Strain  paths  and  stress  paths  at  R  =  2R  cylindrical  wave 
propagation  in  Mixed  Company  sandstone.  A  radial  stress 
given  by  o  =  pQ  exp(-at),  with  (1/a)  =  1  msec  and  p  =  1 
kbar,  is  applied  at  R  =  1  m.  0 


(XI  O'2  ) 


—  ji  »uu  u  •  ju  >uu  — J»w»U 

#1  EPS-T.  »2  CP/K 


Figure  2a.  Strain  paths  and  stress  paths  at  R  =  1.5R0  for  spherical  wave 
propagation  in  Mixed  Company  sandstone.  A  radial  stress 
given  by  a  =  p  exp(-at),  with  1/a  =  1  msec  and  pn  =  10  kbar, 
is  appliedat  R„  =  1  m.  0 


*1  EPS-T, 


#2  CP/K 


CXI  O-3  ) 


Figure  2b.  Same  as  2a,  but  with  R  =  2R  .  Note  changes  in  vertical  and 
horizontal  scales.  ° 


#1  EPS-fl.  #2  L  /  M  U  C  X 1  0"3  ) 

.00  1.50  3.00  4.50  6.00  7.50  9.00  10.50 


Figure  2e.  Same  as  2a,  but  with  R  =  5R  .  Note  changes  in  vertical  and 
horizontal  scales.  0 


R  =  1.5  RQ  R  =  1.5  RQ  R  =  1.5  RQ 

1/a  =  0.1  msec  1/a  -  1-0  msec  l/a  =  iq  msec 


Figure  3.  Strain  paths  and  stress  paths  at  various  positions  for  cylindrical 
wave  propagation  in  Mixed  Company  sandstone.  A  radial  stress 
given  by  o  =  p0  exp(-at),  with  p0  =  10  kbar  and  various  values 
of  1/a,  israpplied  at  R  =  1m.  Note  changes  in  the  vertical  and 
horizontal  scales  in  each  graph. 
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Figure  3.  Continued 


STATIC  EXPERIMENTAL  SIMULATION  OF  LOAD-UNLOAD  PATHS 


Stress  and  strain  paths  were  determined  experimentally  in  the  (L,  pc) 
and  (e^,  e  )  planes  using  the  results  suggested  by  various  one-dimensional 
finite-difference  solutions  given  previously.  A  detailed  discussion  of 
experimental  techniques  used  in  these  tests  is  presented  in  Appendix  II. 

The  stress  and  strain  paths  considered  here  correspond  approximately  to 
those  given  in  Figure  3  for  R  =  3Rq  and  three  separate  decay  constants 
(1/a  =  0.1  msec,  1.0  msec  and  10  msec).  Figures  4a,  4b  and  4c  show  the 
three  characteristic  strain  paths  generated  from  the  numerical  solution  and 
the  strain  paths  to  be  followed  in  the  static  laboratory  tests.  The  percent 
strains  indicated  here  are  used  only  to  indicate  orders  of  magnitude  and  are 
not  the  actual  values  achieved  during  testing.  No  attempt  was  made  to  follow 
the  numerically  determined  strain  paths  exactly,  they  were  used  simply  to 
indicate  the  qualitative  nature  of  load-unload  paths  in  the  vicinity  of 
buried  explosions.  Figure  4a  shows  the  calculated  and  experimental  paths 
corresponding  to  a  decay  time  of  1/a  =  0.1  msec;  this  consists  essentially 
of  uniaxial-strain  loading  and  constant-axial-strain  unloading  followed  by 
uniaxial-strain  unloading.  Figure  4b  shows  the  theoretical  path  corresponding 
to  1/a  =  10  msec  in  comparison  to  the  experimentally  followed  path.  The 
experimental  strain  path  to  be  used  consists  of  a  uniaxial-strain  loading 
and  a  constant-axial-strain  unloading.  Finally  Figure  4c  shows  the  theoretical 
path  for  1/a  =  10  msec  in  comparison  to  the  experimentally  followed  path. 

The  experimental  strain  path  to  be  used  consists  of  uniaxial-strain  loading 
and  constant-volume-strain  unloading.  Kayenta  sandstone  from  the  Mixed 
Company  site  was  the  material  tested  in  this  investigation. 
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uniaxial-strain  unloading.  The  unloadinq.  The  percent  strains 

percent  strains  indicated  here  are  indicated  here  are  used  to  indicate 

used  only  to  indicate  orders  of  orders  of  magnitude  and  are  not  the 

magnitude  and  are  not  the  actual  actual  values  achieved  during  testi 

values  achieved  during  testing. 
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gure  4c.  Comparison  of  the  theoretical 

(calculated)  strain  path  to  the 
experimental  strain  path  to  be  followed 
during  testing  of  path  III  (1/a  = 

10  msec).  The  experimental  path 
shows  a  uniaxial-strain  loading 
with  a  constant-volume-strain 
unloading.  The  percent  strains 
indicated  here  are  used  only  to 
indicate  orders  of  magnitude  and 
are  not  the  actual  values  achieved 
during  testing. 
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TEST  RESULTS 


The  three  strain  paths,  I,  II  and  III,  used  in  testing  the  Kayenta 
sandstone  are  shown  in  Figures  5a,  6 a  and  7a,  respectively.  Since  all 
loading  was  conducted  under  uniaxial -strain  conditions,  a  composite  loading  curve 
is  shown  for  each  path  type.  Individual  unloading  curves  are  shown  for 
each  test,  departing  from  the  composite  loading  curve  at  their  respective  maximum 
strains.  The  stress  paths  generated  from  the  three  strain  paths  are  shown 
in  Figures  5b,  6b  and  7b.  Composite  loading  curves  are  shown  along  with 
individual  unloading  curves.  Included  in  each  stress  path  figure  is  the 
triaxial  failure  envelope  generated  from  this  material.  Tables  I,  II  and 
III  give  computer  listings  for  each  test.  Table  Column  1  gives  the  data 
point  while  columns  2  through  8  give  confining  pressure  (pc)  in  kilobars, 
axial  load  (a  - p  )  in  kilobars,  axial  strain  (el  in  percent,  the  two 

a  C  a 

transverse  strains  (e.  and  e*  )  in  percent,  volume  strain  (ea  +  e.  +  ) 

1*1  j*2  d  l2 

in  percent  and  mean  stress  [1/3(ct  +2p  1]  in  kilobars.  All  plots  were 

a  C 

constructed  from  these  tables. 


LATERAL  STRAIN,  % 

Figure  5a.  Strain  path  followed  during  uniaxial-strain 
loading  and  constant-axial  and  uniaxial- 
strain  unloading. 


STRESS  PATH  FOLLOWED  DURING 
STRAIN  PATH  I  TESTING 
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’  iutp  St.  Stress  path  followed  during  uniaxial- 
strain  loading  and  constant-axial 
and  uniaxial-strain  unloading.  The 
resulting  stress  path  is  a  composite 
of  four  tests. 
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Figure  6a. 


Strain  path  followed  during  uniaxial 
strain  loading  and  constant-axial- 
strain  unloading.  The  resulting 
stress  path  is  a  composite  of  four 
tests . 
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Figure  6b. 


Stress  path  followed  during  uniaxial- 
strain  loading  and  constant-axial- 
strain  unloading.  The  resulting 
stress  path  is  a  composite  of  four 
tests . 
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Figure  7a.  Strain  path  followed  during 
uniaxial -strain  loading  and 
constant-volume-strain  un¬ 
loading.  The  resulting  stress 
path  is  a  composite  of  three 
tests. 
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Figure  7b.  Stress  path  followed  during 
uniaxial- strain-loading  and 
constant-volume-s train  unload¬ 
ing.  The  resulting  stress 
path  is  a  composite  of  three 
tests. 
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V 


*  Axial  strain  rezeroed  for  constant-axial -strain  unloading 

*  Lateral  strains  rezeroed  for  uniaxial -strain  unloading. 
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Axial  strains  rezeroed  for  constant-axial-strain  unloading. 
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Axial  strain  rezeroed  for  constant-axial -strain  unloading 
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Shows  only  the  uniaxial -strain  loading. 


DISCUSSION  AND  CONCLUSIONS 


An  initial  observation  of  the  experimentally  observed  stresses 
indicates  that  there  is  little  difference  between  the  load-unload  paths 
for  strain  path  types  II  and  III.  Figures  6b  and  7b  suggest  that  if  the 
yield  condition  is  reached  during  uniaxial-strain  loading  with  the  stress 
paths  following  along  the  yield  surface,  then  the  unloading  stress  paths 
are  similar  in  direction  and  magnitude  for  either  constant-axial -strain  or 
constant-volume-strain  unloading.  The  numerial  analysis  solutions  agree 
with  the  above  observation  in  that  regardless  of  the  strain  path,  the 
stress  path  would  follow  along  the  yield  surface  during  unloading  (provided 
that  yield  was  reached  during  uniaxial -strain  loading).  All  of  the 
experimentally  observed  stress  paths  show  the  unloading  curve  to  go  initially 
above  and  then  cross  through  and  go  below  the  loading  curve.  The  experimental 
unloading  curves  did  not  remain  on  or  intersect  (as  in  the  case  of  strain 
path  2)  the  yield  surface  as  illustrated  by  the  numerical  analysis. 

Such  variations  in  unloading  material  behavior  may  be  modeled  by 
including  additional  phenomena  into  the  constitutive  equations.  Phenomena 
to  be  included  in  the  equations  would  be  permanent  volume  compaction  and 
work-hardening  of  the  shear  failure  envelope.  The  former  effect  will 
mainly  influence  the  strain  paths  and  the  latter  will  change  the  stress 
paths,  particularly  in  the  unloading  portion.  It  was  experimentally 
determined  that  the  material  behaved  nonlinearly  during  initial  loading 
as  compared  to  the  linear  model  used  in  the  numerical  analysis.  Such 
nonlinearities  may  be  also  handled  by  the  aforementioned  considerations. 

The  observation  that  the  unloading  path  lies  below  the  loading  path  in 
stress  space  may  be  related  to  fracture  and  the  resulting  loss  of  cohesion, 
rather  than  ductile  plastic  flow,  as  asumed  in  the  calculations. 
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Inclusion  of  pore  pressure  effects  into  the  model  would  be  of  interest 
in  future  work.  Both  the  calculations  and  laboratory  strain-path  tests 
should  be  performed  under  various  saturation  conditions.  Much  of  the  previous 
theoretical  work,  including  the  finite-difference  computer  code,  already 
contains  this  capability;  it  has  just  not  been  exercised  yet.  Also  of  future 
interest  would  be  some  theoretical  results  for  two-dimensional  dynamic  loading 
situations,  expressed  in  terms  of  e  .  £. ,  l  and  p  .  This  could  be  done 

a  L  C 

by  calculating  the  following  invariants  as  functions  of  time  at  a  particular 
material  element: 


T  ( t ) 

=  (1/6)[(oji 

-a22)2  +  (a22-o33)2  +  (a33-0ll)2]  4 

)V2 

Oi22  + 

0 1 32  4 

°232  | 

04) 

P(t) 

=  (oil 

+  a22 

+  cr  3  3  )/3  , 

05) 

*v(t) 

=  EH 

4  e22 

4  e33 

06) 

cd(t)  =  j  (l/6)[(cn-e22)2  +  (e22-e33)2  +  (e33-Ell)2]  + 

0/2 

M22  +  c13"  +  C232  .  (17) 


The  desired  quantities  used  for  comparison  with  laboratory  tests  are  then 
obtained  from  Eqs.  (9)  -  (12). 

The  results  presented  here  have  shown  that 

(1)  We  can  define  strain  paths  for  static  testing  of  rock  (and  soil) 
samples  that  are  more  representative  of  actual  field  situations  that 
those  commonly  used  heretofore  in  constitutive  modeling,  and  that 

(2)  It  is  possible  to  reproduce  these  paths  in  laboratory  tests. 
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APPENDIX  I 


General  Re  I ationships  and  Finite-Difference  Calculations 

The  equation  for  momentum  conservation  in  Eulerian  coordinates  is 
given  by 


v 


+ 


(18) 


where  ,  is  the  material  density,  v  is  the  radial  particle  velocity,  o 
and  o  are  the  radial  and  tangential  stress  components,  and  g  is  1  (for 
plane  flow),  2  (for  cylindrical  flow)  or  3  (for  spherical  flow).  A  dot 
over  a  variable  indicates  time  differentiation  at  a  fixed  material  element 
and  r  is  the  Eulerian  spatial  coordinate.  It  is  inconvenient  to  deal  with 
Eulc?  an  coordinates,  hence  we  choose  to  express  Eq.  (18)  in  terms  of 
Lagrangian  coordinates  representing  the  initial  configuration.  We  define 
R  as  the  initial  radial  coordinate  of  a  material  element  whose  current 
radial  location  is  at  r.  Radial  and  transverse  stress  components  in  the 
initial  configuration  (Lagrangian)  are  denoted  aD  and  a..  If  the  initial 
density  is  given  by  pQ,  then  mass  conservation  requires  that 


DQRg'1  dR  =  pr9’1  dr 


09) 


If  the  forces  on  a  material  element  are  to  be  the  same  in  the  two 
representations,  then 


(20) 
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Now  write  Eq.  (18)  as 


-pr9'^  dr  ^  =  d(r9"^o  )  -  o  dr9  "* 

r  w 


(22) 


keeping  in  mind  that  the  differentials  on  the  right-hand  side  are  taken 
at  constant  time.  Substitution  of  Eqs.  (19)  -  (21)  into  Eq.  (22)  then  gives 


-PoR9'1  dR  v  =  d(R9'1aR)  -  a^R9"1 


or 


3a  c 

-V  =  w 


°R  -  °0 


(23) 

(24) 


in  Lagrangrian  coordinates. 

In  order  to  use  Eq.  (24)  in  a  finite-difference  solution,  an  artificial 
viscous  stress  q  is  included.  The  following  equations,  with  the  addition 
of  a  constitutive  law,  then  form  the  basis  of  the  numerical  calculations: 


p 


3or 

3  R 


(g 


(25) 


q 


3  V|  2 


p0a2  (ar)2  m 


3  V 

3  R 


S  0 


=  0 


3V 

3R 


>  0 


3  V 

3  R 


(26) 


(27) 


where  A  is  nondimensional  constant  on  the  order  of  unity,  AR  is  the  spatial 

increment  in  the  finite-difference  solution,  and  eR  and  e'Q  are  the  radial 

and  tangential  strain  rates  in  the  initial  configuration.  A  straight-forward 

centered  difference  scheme  is  used  and  Eqs.  (25)  -  (27)  are  written  in 
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finite-difference  form  as 


i+%  i 
vj  '  vj 

At 


(On)’  -  ( 0  p ) 1 

R  i+h  R  i-H 
aR 


(od)1'  +  (Op)’  -  (o  )]  -  (a  )\ 

(g-1)  j+h _ j-h _  j+h _ j-*s 

2Ri 


i -H  i -h 
*  qi-* 


AR 


(28) 


Ur)1** 

j+h 


v1+2  -  v1:2 

_ j±i 

AR 


(29) 


i+%  .  i  +!j 

(e  )i+%  =_!i _ _!m 


(30) 


The  stress  rates  (aR  and  aQ)  are  obtained  from  eR  and  eQ,  and  therefore 
the  stresses  and  strains  are  calculated  from 


(31) 


where  X  represents  oR,  aQ,  eR  and  e0. 

The  constitutive  model  used  here  is  expressed  in  terms  of  the  principal 
stress  and  strain  components  and  (i  =1,2  and  3)  with  the  following 
identification: 
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g 


1  (Plane  Flow) : 


£^  —  -3v/3R  9  6^2  =  £  3  ”  0 

0 1  =  cr^»  O2  =  03  = 

g  =  2  (Cylindrical  Flow) 

£1  =  -9v/aR  ,  £2  =  -v/R  ,  £3=0 

°1  =  °R  >  a2  =  aQ  »  a3  =  al 

g  =  3  (Spherical  Flow): 

£l  —  -8v/ 3R  ,  £2  —  £3  =  -v/R 

aI  =  aR  *  °Z  =  °3  =  °q 

Let  us  define  the  volume  strain  £y,  the  mean  stress  p,  the  stress 


deviators  s^ 

and  the  second  invariant  of  the  stress 

tenser  according  to 

ev 

=  ei  +  £2  +  e3 

(32) 

P 

=  ( a  x  +  o2  +  a3)/3  , 

(33) 

si 

Cl 

I 

•r— 

D 

II 

(34) 

J2 

=  (s2  +  s2  +  s2)/2 

(35) 

The  elastic-plastic  constitutive  relation  used  here  is  then  defined 
according  to  the  following  equations: 


40 


m 


P  =  PUV)  »  (36) 

s,  =  2 u(e.  -  c  /3)  -  2UC—  •  (37) 

1  TV  /jj 

The  variable  £  is  determined  by  the  condition  that  the  stress  state  must 
remain  on  the  failure  surface,  defined  by 

^37  =  f(P)  .  (38) 

when  a  material  element  is  undergoing  plastic  deformation. 

From  Eq.  (35)  we  find  that 

2/37  ^2  ~  sjsi  (Summation)  (39) 

and 

ST2  =  (u//J7)  s-e.  -  2uc  =  f'(p)p  .  (40) 

Therefore,  the  variable  £  in  Eq.  (37)  is  given  by 

2p£  =  (p//37)  s1ei  -  f'(p)p  ,  (41) 

or,  in  terms  of  a..  and  p,  as 

2p£  =  (u/v/7)(o.e.  -  pev)  -  f'(p)p  .  (42) 

If  it  is  desired  to  include  effects  of  fluid  saturation  defined  by 
nonzero  pore  pressure  Pp,  a..  is  replaced  by  the  effective  stress  components 
<0.>  =  a.  -  nPp  (0  <  n  <  1)  in  the  elasticity  relationship  and  by  a..*  h  -  Pp 
in  the  failure  surface  ralationship: 

<p>  =  p  -  nPp  =  p(ey)  ,  (43) 


(44) 


'  *  ^  i 

<s^  =  Si  =  2p ( -  ev/3)  -  2V^-— 

/Jo 


2y£  =  (y//37)(a.e.j  -  pev)  -  f 1  (p*)  ( 1 -m)|> 


where 


dp 


The  function  f(p)  is  taken  to  be  of  the  form 


f(p)  =  SQ  +  AS(1  -  e'p/a) 


(45) 


(46) 


(47) 


Analytical  Determination  of  Elastic  Stress  and  Strain  Paths  for  a 
Spherical  Explos ion 

If  u(r,t)  is  the  radial  displacement,  the  spherical  wave  equation 
for  purely  elastic  deformation  can  be  written  as 

32u/3t2  =  c2[32u/3r2  +  (2/r)  3u/3r  -  (2/r2)u]  ,  (48) 

where  r  is  the  radial  coordinate,  t  is  the  time  and  c  is  the  longitudinal 
elastic  wave  speed.  This  expression  takes  a  simpler  form  if  it  is  written 
in  terms  of  a  displacement  potential  f  such  that 

u(r,t)  =  c23/3r  (Y/r)  .  (49) 

In  this  case 

32H'/3t2  =  C2  324l/3r2  ,  (50) 

whose  solution  for  outgoing  waves  is  given  by  the  familiar  expression 

»  =  f(t - ~)  .  (51) 
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The  displacement,  strain  components  and  stress  components  can  be  expressed 
in  terms  of  y  and  its  derivatives  y'  and  y"  according  to 


u(r,t) 

=  -(c/r)y‘  -  (c/r)2  y  , 

(52) 

~ea  = 

3u/3r  =  (l/r)y“  +  (2c/r2)y*  +  (2c2/r3)y  , 

(53) 

'et 

u/r  =  -(c/r2)y‘  -  (c2/r3)y  , 

(54) 

_0a 

(1/r)  [(A+2M)y"  +  {4yc/r)y'  +  (4yC2/r2)y]  , 

(55) 

-°t  = 

(1/r)  [Ay"  -  (2uc/r)y‘  -  (2yc2/r2)y] 

(56) 

where  A  and  y  are  the  Lame  constants.  The  sign  convention  used  throughout 
this  work  is  that  stresses  and  strains  are  positive  in  compression.  For 
a  pressure  history  at  r  =  rQ  given  by 


or(rQ,t)  =  0 


t  <  0 


°r^ro,t^  =  poe 


t  ^  0 


The  function  y  must  satisfy  the  following  ordinary  differential  equation: 

(A+2y)f"(t)  +  (4yC/ro)4',(t)  +  (4Pc2/ro2)y(t)  =  (58) 


•r  p„e 
oro 


subject  to  the  conditions,  from  Eqs.  (52)  and  (58),  that  jumps  in  y 
and  y*  at  t  =  0  obey  the  following  relationships: 


(A  +  2p)  O']  +  (4uc/r0)  [y]  =  0 


O']  +  (c/ro)  [y]  =  0 
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where  [  ]  indicates  the  jump  in  the  function,  i.e.,  [f]  =  f(0+)  -  f(0"). 
Equations  (59)  thus  require  that  f  and  each  be  continuous  at  t  =  0  as 
long  as  a  f  2p.  Hence,  a  solution  to  Eq.  (53)  can  be  written  as 


where 


y(t)  =  e~S2t  (M  cos  b j t  +  N  sin  6}t)  +  Yoe~at 


(60) 


M  =  -v_  = 


r  P 
o*o 


0  a2(x+2y)  -  4uCa/rQ  +  4pC2/rQ2 


(61) 


ar  (a+2m)  -  2uc 

N  =  -■  - - v 

2c  ^u(a+u) 


(62) 


a  =  2C 

01  rQ  (X+2UJ 


(63) 


r  -  2uc 

02  _  rQ  (I+2„) 


(64) 


In  the  case  of  an  elastic  fluid  u  =  0  and  the  displacement  potential 
and  its  first  two  derivatives  become 


*  .  ga  <1  -  -  at) 


(65) 


r  p 

(e-°t  .  ,) 

Act 


(66) 


V '  -  - 


Vo  e-at 


(67) 


If  «  =  0  (i.e.,  the  cavity  pressure  remains  constant  at  p  )  in  the  case  of 


a  fluid,  the  displacement  potential  and  its  first  two  derivatives  become 


v  =  '  ~2ir  t2  *  (68) 

f*  *  -  ~~  t  •  (69) 

r  p 

*"  =  -  ~~  .  (70) 

In  the  special  case  of  spherical  wave  propagation  we  can  make  the 
identification  that  L  =  a,  -  a.  and  p„  =  a.,  in  which  case  the  stress  and 

a  l  C  L 


strain  paths  can  be  written  parametrically  as 

L  =  -(2M/r)0"  +  Oc/r)'?'  +  (3c2/r2)y]  ,  (71) 

Pc  =  -0/r)[Ar  -  ( 2uc/r )?’  -  (2u c2/r2)?]  ,  (72) 

ea  =  -(l/r)[»"  +  (2c/r)v'  +  (2c2/r2h]  ,  (73) 

u 

et  =  (c/r2)  [V  +  (c/r)4']  .  (74) 


Equations  (71)  to  (74)  in  the  case  of  spherical  elastic  waves  are  the 
analytical  counterparts  of  Eqs.  (9)  to  (12)  for  numerical  solutions. 
Comparison  of  strain  and  stress  paths  calculated  by  the  two  methods  is 
shown  in  Figure  8  for  1/a  =  1  msec,  R/RQ  =  3,  K  =  95  kbar,  c  =  3  km/sec, 
and  pQ  =  2.0  gm/cm3.  It  can  be  seen  that  the  numerical  solution  gives  a 
good  approximation  of  the  strain  and  stress  paths  except  for  the  peak 
values  associated  with  the  main  compressive  fronts.  This  is  a  result  of 
the  viscous  stresses  that  are  included  in  the  finite-difference  solution  to 
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damp  out  numerical  oscillations,  and  has  no  significance  with  regard  to 
the  conclusions  reached  in  this  report. 


APPENDIX  II 


EXPERIMENTAL  TECHNIQUE 


Specimen  Preparation 

Specimens  were  prepared  from  Kayenta  sandstone.  Mixed  Company  Site. 
Cylindrical  samples  3.81  centimeters  long  by  1.91  centimeters  diameter 
were  used  thus  maintaining  a  length  to  diameter  ratio  of  2  to  1 .  Specimen 
ends  were  ground  parallel  to  within  ±  .001  centimeters .  Specimens  were 
air  dried  with  weight,  length  and  diameters  being  recorded  for  each  sample 
for  use  in  determining  sample  density  and  strains.  Samples  were  prepared 
for  testing  by  first  wrapping  them  in  urethane  plastic  (.025  cm  thick) 
with  hardened  steel  endcaps  attached  at  each  end  using  stainless  steel  lock 
wire. 

Stress  and  Strain  Determination 

Stress  and  strain  transducers  were  placed  within  the  pressure  vessel. 
Confining  pressure  was  measured  using  a  calibrated  350-ohm  manganin  pressure 
sensitive  coil  accurate  to  ±  .003  kbars.  Jacketed  samples  were  placed 
and  centered  on  the  load  cell  when  in  the  pressure  vessel.  The  load  cell 
was  accurate  to  +  .005  kbars.  Axial  and  lateral  strain  transducers  were 
of  the  cantilever  type  using  strain  gauges  in  a  Wheatstone  bridge  configuration 
to  obtain  voltage  output.  The  axial  cantilevers  measured  total  axial 
displacement  and  were  calibrated  to  be  accurate  to  +  .003  percent  strain. 
Lateral  strain  cantilevers  were  positioned  at  mid-sample  and  sampled  strains 
at  90  degree  intervals.  Diametrically  opposed  arms  were  calibrated  for  lateral 
strain.  The  lateral  strains  were  averages  with  a  resulting  accuracy  of  ±  .006 
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percent.  Figure  9  shows  a  schematic  of  the  transducers  when  inside  the 
pressure  vessel.  Further  discussion  on  transducer  design  may  be  obtained 
in  Terra  Tek  report  TR  75-29. 

Testing  Procedures 

Seven  samples  were  first  tested  triaxially  to  failure  to  generate  the 
triaxial  failure  envelope  for  the  material  while  eight  samples  were  tested 
following  the  three  strain  paths.  Triaxial  testing  commenced  by  first 
hydrostatically  loading  the  samples  to  the  desired  confining  pressure  with 
subsequent  axial  loading  to  failure,  stresses  and  strains  being  recorded 

-4  -  1 

during  all  phases  of  loading.  A  strain  rate  of  about  10  sec  was  used 
during  loading. 

Uniaxial-strain  loading  was  used  when  following  a  specified  strain 
path.  Axial  load  and  confining  pressure  were  applied  such  that  zero 
lateral  strain  was  maintained.  When  following  strain  path  I,  II  or  III 
during  unloading,  i.e.,  constant-axial -strain  and  uniaxial -strain  unloading, 
constant  axial  strain  unloading  and  constant  volume  strain  unloading, 
respectively,  the  confining  pressure  and  axial  load  were  adjusted  to 
maintain  the  desired  strain  state. 

Data  Acquisition  and  Analysis 

Both  x-y  recorders  and  a  POP  Lab  11  computer  were  used  for  data 
acquisition.  The  x-y  recorders  were  used  primarily  for  instantaneous 
feedback  during  testing  while  the  PDP  Lab  11  computer  data  was  used  for 
analysis  of  pressure  effects,  endcap  effects  and  generation  of  stress  and 
strain  load-unload  curves.  Tables  I,  II  and  III  presented  in  the  text 
are  a  result  of  the  computer  analysis. 
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AXIAL  STRAIN  CANTILEVERS 
LATERAL  STRAIN  CANTILEVERS 


-  CONFINING  PRESSURE  PORT 


Figure  9.  Pressure  vessel  schematic  showing  the  sample  and 
stress  and  strain  transducers. 
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DIFFERENTIAL  STRESS  (  CTj  -CT3  )jKbar 


Figure  9a.  Stress  path  followed  during  strain  path  III  testing. 
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CONFINING  PRESSURE  (03 ),  kbar 
Figure  9c.  Stress  path  followed  during  strain  path  II  testing. 
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DIFFERENTIAL  STESS  (CT-  <r,),  kbar 
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Figure  9f.  Strain  path  followed  during  path  III  testing 
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